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Abstract 

■ We provide a new natural interpretation of the Lax representation 

, for an integrable system; that is, the spectral problem is the linearized 

form of a Miura transformation between the original system and a 
modified version of it. On the basis of this interpretation, we formu- 
late a systematic method of identifying modified integrable systems 
I that can be mapped to a given integrable system by Miura transfor- 

QQ ■ mations. Thus, this method can be used to generate new integrable 

. systems from known systems through inverse Miura maps; it can be 

applied to both continuous and discrete systems in 1 + 1 dimensions 
as well as in 2 + 1 dimensions. The effectiveness of the method is illus- 
trated using examples such as the nonlinear Schrodinger (NLS) sys- 
tem, the Zakharov-Ito system (two-component KdV), the three- wave 
interaction system, the Yajima-Oikawa system, the Ablowitz-Ladik 
lattice (integrable space-discrete NLS), and two (2 + l)-dimensional 
^ . NLS systems. 
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1 Introduction 



For an integrable system that is essentially nonlinear and thus is not lineariz- 
able by a change of variables, the Lax (or zero-curvature) representation [1] 
occupies a central position in its integrability properties; giving the Lax rep- 
resentation is even considered as proof of integrability. The first example of 
a Lax representation was found in the late 60s for the Korteweg-de Vries 
(KdV) equation [2]; it was originally derived through linearizing a Riccati- 
type transformation from the modified KdV (mKdV) equation (or, more 
generally, its one-parameter generalization called the Gardner equation [3]) 
to the KdV equation [HIS]. Such a nontrivial and non-ultralocal transfor- 
mation of dependent variables from one integrable system to another is now 
called a Miura map P,!^. Subsequently, in the early 70s, the Lax represen- 
tation in 2 X 2 matrix form was found for the nonlinear Schrodinger (NLS) 
equation |6], without any reference to the relevant Miura map. Since then, 
the number of integrable systems admitting the Lax representation has been 
increasing rapidly, including multicomponent systems, higher dimensional 
systems [3-|9], and discrete systems [lO]. However, as the Lax representa- 
tions have been extended in various directions, the primary role played by 
the original Miura map in the KdV case has gotten lost in oblivion. 

The main theme of this paper is to demonstrate, using an abundance 
of specific examples, that the Miura maps are by no means less important 
than the Lax representations. In fact, a Miura map can generally produce 
the relevant Lax representation and vice versa. Thus, they are different 
facets of the same property and, in a sense, equivalent. This is not merely 
a conceptual integration of the Miura maps and the Lax representations; 
rather, we will use it for a more practical and attractive purpose, that is, to 
generate new integrable systems from known systems in a systematic manner. 
Actually, no such methods have ever been proposed and applied successfully 
to a broad spectrum of examples. The basic idea of our new method can be 
best described using the nonreduced NLS system [11] 



as an illustrative example. The subscripts t and x denote the partial differ- 
entiation with respect to these variables. Note that the NLS system (11. ip 
is integrable for matrix- valued dependent variables [12]; in the general case, 
Q and R are li x I2 and I2 x li matrices, respectively. In this paper, the 
symbol O on the right-hand side of the equations implies that the dependent 
variables can take their values in matrices. The Lax representation for the 



iQt + Q.. - 2QRQ = O, 
iRt — Rxx + 2RQR = O, 
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matrix NLS system (II. ip is given by [T^fTl] 
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(1.2a) 
(1.2b) 



Here, C, is the spectral parameter, which is an arbitrary constant independent 
of X and t, and Ji and I2 are the li x li and I2 x I2 unit matrices, respectively. 
We consider an (/i + ^2) x ^1 matrix- valued solution for the pair of linear 
equations (11.21) such that is an /i x li invertible matrix. Then, in terms 
of the I2 X /i matrix P := ^s^FS (Q can be rewritten as a pair of matrix 
Riccati equations (see P^5HT7] for the scalar case), 

P^ = R + 2iCP - PQP, (1.3a) 
Pt = 2CR - iR^ + 4iC^P + iRQP + iPQR - 2CPQP - iPQ^P. (1.3b) 

Using fll.3ap . we can express R in terms of P and Q as — 2i^P + P^ + PQP. 
Thus, (11. lap and fll.3bp now comprise a closed system for Q and P, i.e., 

iQt + Qxx + ^KQPQ - '2QPxQ - 2QPQPQ = O, (1.4a) 
iPt - P,, - AiCPQP - 2PQ,P + 2PQPQP = O. (1.4b) 

This is intrinsically a derivative NLS system, which was investigated by 
Ablowitz et al. [TS] and Gerdjikov and Ivanov [T^] in the case of scalar de- 
pendent variables. The matrix generalization (II. 4p was studied in [SUHSl] . 
Note that the free parameter C, is nonessential as long as we consider the 
Gerdjikov-Ivanov system (ll.4p separately as an isolated system; indeed, it 
can be set equal to zero using a Galilean transformation (cf. [25l|26]). It 
is straightforward to confirm that if the pair {Q,P) satisfies (11.41) . then the 
pair (Q, R) with R := -2iCP + Px + PQP indeed satisfies (dH]); this fact is 
already known in the case C = (see [TM7rH2U] for the scalar case and [20tl25] 
for the matrix case). Thus, starting with the NLS system (11.11) and its Lax 
representation (11.21) . we obtain a Miura map from the Gerdjikov-Ivanov sys- 
tem (II. 4p to the NLS system (II. ip . without any prior knowledge about (11.40 . 
The spectral problem (ll.2ap turns out to be the linearized form of the defining 
relation (ll.3ap for the Miura map (cf. [501151] ): the entire Lax representation 
(11.21) not only yields the NLS system (11.11) as the compatibility condition but 
also enables us to apply an inverse Miura map to (II. ip . Somewhat similar 
computations were given in [321. 133] . but they are not as straightforward and 
comprehensible as ours. The "nonessential" parameter ( plays a crucial role 
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when studying the original system f ll.ip and the modified system (ll.4p in a 
unified way. For instance, the simplest conservation law for f ll.4p . 



i-tr(gP) + — tr [Q,P - QP, - (QPf] = 0, (1.5) 

works as a generating function of an infinite set of conservation laws for 
(11. ip . Indeed, following the same routine as in the case of the original Miura 
map [SlIS]' can formally express P using (ll.3ap in a power series of {2iQ~^, 

„ R Rr Rrr — RQR ,^ „, 

P = — \ . (1.6) 

2iC (2iC)2 (2iC)3 ^ ^ ^ 

Thus, substituting (II. 6p into (II. 5p and equating all the coefficients of differ- 
ent powers of (2iC)~^ with zero, we obtain an infinite set of conservation laws 
for the NLS system (11.10 : this uncovers the hidden essential meaning of the 
method devised by Wadati et al. [17] (also see [31]) for constructing conserva- 
tion laws recursively. Then, using (ll.Sap . namely, R = —2i(P + P^ + PQP, 
we can also obtain an infinite set of conservation laws for the Gerdjikov- 
Ivanov system (II. 4p . Moreover, using the inverse scattering method based 
on the Lax representation (II. 2p . it is possible to derive the solution formulas 
for the NLS system (II. ip and for the Gerdjikov-Ivanov system (11.41) concur- 
rently. 

Our method described above actually applies to a surprisingly wide spec- 
trum of integrable systems and is not limited to systems of two coupled 
partial differential equations (PDEs) in 1 + 1 space-time dimensions, such as 
(11.11) . Broadly speaking, the method is applicable to every integrable system 
with a Lax representation, the spatial part of which is ultralocal in the depen- 
dent variables; then, the method successfully identifies modified integrable 
systems embedded in the original system via Miura maps. Thus, it can be 
applied to most of the known integrable systems, including multicomponent 
PDEs, (2 + l)-dimensional PDEs, differential-difference equations, and par- 
tial difference equations. The Lax representation for an integrable system can 
be determined only up to gauge transformations; we must fix the gauge ap- 
propriately to write the spectral problem in an ultralocal form. For instance, 
the Gerdjikov-Ivanov system (ll.4p with C = Ci allows the nonstandard Lax 
representation, namely, (ll.2p with R = —2i(iP + Px + PQP] this is not a 
convenient form for applying our method, because the spectral problem is 
not ultralocal in P. Thus, before applying the method, we consider a suit- 
able gauge transformation, e.g., $i := \E'i and $2 '■= ^2 ~ -P^i, so that the 
spectral problem assumes an ultralocal form with respect to both Q and P. 
However, such a procedure is not always possible for a nonstandard spectral 
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problem and our method can, in general, be applied to a given integrable sys- 
tem with only a finite number of iterations. That is, it terminates at some 
stage, resulting in a finite chain of j-th modified systems (j = 0, 1, . . . , jmax)- 
The remainder of this paper is organized as follows. In section 2, we de- 
scribe the method in a general setting. First, we reinterpret the spatial part 
of a given Lax representation as the defining relations of a Miura map; sub- 
sequently, the temporal part of the Lax representation is used to emboss an a 
priori unknown modified system related to the original system by the Miura 
map. For an integrable system in 1 -|- 1 dimensions, the spatial part of the Lax 
representation is genuinely a spectral problem, that is, it involves intrinsically 
the spectral parameter, which is an arbitrary constant. Thus, the relevant 
Miura map and modified system always contain an additional free parameter 
that the original system does not have. Using this parameter as an expansion 
parameter, the simplest conservation laws for the modified system can gen- 
erate infinitely many conservation laws for the original system, which in turn 
provide infinitely many conservation laws for the modified system. Sections 3 
and 4 are devoted to specific examples of such (1 -|- l)-dimensional integrable 
systems. In section 3, we apply the method to continuous systems, namely, 
PDEs such as the Zakharov-Ito system p!3|[35] , which is a two-component ex- 
tension of the KdV equation, the Jaulent-Miodek system |36] , the three- wave 
interaction system [37], and the Yajima-Oikawa system [38]. In section 4, 
we investigate space-discrete systems; the Toda lattice in Flaschka-Manakov 
coordinates [55H1T| . the Belov-Chaltikian lattice [12], the relativistic Toda 
lattice [13], and the Ablowitz-Ladik lattice ^1] are discussed as illustrative 
examples. In section 5, we consider integrable PDEs in 2 -f 1 dimensions. 
Their Lax representations can be divided into two distinct groups according 
to the dimensionality of the spatial Lax operators [13]. In the first group, the 
spatial part of a Lax representation is a one-dimensional spectral problem as 
in the (1 -|- l)-dimensional case [151 - H7] : only the associated time evolution ex- 
hibits a (2 -|- l)-dimensional nature. In the second group, the spatial part of a 
Lax representation is genuinely a two-dimensional problem; in contrast to the 
(1 + l)-dimensional case, it does not contain any essential spectral parame- 
ter. Our method can be applied to both groups, but in the latter case, the 
resulting Miura map and modified integrable system contain no additional 
parameter. We take one representative example from each group, namely, 
the Calogero-Degasperis system [16] and the Davey-Stewartson system [18], 
to illustrate the method; they are both (2 + l)-dimensional generalizations 
of the NLS system fll.ip . Actually, our method is also applicable to discrete 
systems in 2 + 1 dimensions. However, the examples that we know appear 
to be too complicated to discuss here. Section 6 is devoted to concluding re- 
marks. In the appendices, we present a natural spin-off version of our method 
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using illustrative examples. In appendix \^ we prove that the new pair of 
dependent variables {'^^^Q, "^2) defined from the Lax representation fll.2p 
for the NLS system (11. II) satisfies a closed system, which is equivalent to a 
derivative NLS system, called the Chen-Lee-Liu system [19]. In appendix [Bl 
and appendix [Cl we obtain analogous results for the discrete case and the 
(2 + l)-dimensional case, respectively. 

2 Method: from Lax representation to Miura 
map 

In this section, we describe the essence of our method using a general Lax 
representation in which the spatial part is a one- dimensional matrix spectral 
problem that is first order in the space variable. It can be easily modi- 
fied/extended to the case where the spatial part of a Lax representation is 
a higher-order scalar problem or a two-dimensional matrix problem; some 
specific examples will be considered in subsequent sections. 

We consider a pair of linear equations for a column-vector function i/j, 

fiV = UiC)^, = \/(C)V- (2.1) 

Here, Tj {j = 1, 2) are linear operators and U{Q and V{() are square matrices 
depending on the spectral parameter (. Note that in the (2 + l)-dimensional 
case, C is not required to be a constant; it is allowed to depend on some 
independent variables, as illustrated in subsection 15.11 The first equation 
in (12.1 p is the spatial part of the Lax representation, wherein the linear 
operator Ti denotes the partial differentiation by x (Tii/j := d^il^) in the 
continuous-space case or the forward shift operator (Ti'j/'n := V'n+i) the 
discrete-space case. Similarly, the second equation in (12. ip is the temporal 
part of the Lax representation; in the (1 -|- l)-dimensional case, T2 denotes the 
partial differentiation by t (T21/? := dftl^) in the continuous-time case or the 
forward shift operator (T2'j/'„m •= V'nm+i) ^^e discrete-time case. In the 
(2 + l)-dimensional case, T2 denotes a suitable linear combination of partial 
differential operators, e.g., T2 = dt + f{C)dy [SDH52]. It is more convenient 
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to write (12. II) explicitly in the component forms 
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where / > 2. The entries of the matrices U and V are classical quantities, 
i.e., Uij and Vij do not involve operators; however, they can take their values 
in submatrices if U and V are broken into blocks. In such a case, we collect 
and align a set of linearly independent column-vector solutions of (12. 2p so 
that some j become invertible matrices if needed (see below). 

The compatibility condition T1T21P = T2T11P for the overdetermined sys- 
tem (12.11) results in some nontrivial relation between U and V; this relation 
is often referred to as the zero-curvature condition. For example, in the con- 
tinuous (1 + l)-dimensional case, we take Ti = dx and T2 = dt to obtain the 
equation [SSj lSl] 

Ut-Vx + UV-VU = o. 

If we specify the ^-dependent matrices U and V appropriately, the zero- 
curvature condition provides a closed nonlinear system for some (^-independent 
quantities in U and V (cf. (II. ip and (II. 2p ). In such a case, (12. ip gives the 
Lax representation for the integrable system. For brevity, we assume that 
the spatial Lax matrix U is ultralocal in the dependent variables; that is, if U 
involves some dependent variable, then U does not involve its derivatives and 
shifts with respect to the independent variables. In this section, we express 
the entire set of functionally independent dynamical variables contained in 
U as {gi,g2, • • • An}] these variables may take their values in matrices. 
Although our method is also applicable to nonevolutionary systems that are 
negative flows of integrable hierarchies, for notational convenience, we dis- 
cuss the simpler case of evolutionary systems. Thus, the time evolution of 
the integrable system considered can be written explicitly as 



dtQi = fi{qi,q2, - ■ ■ ,qN), i = l,2,...,N 
in the continuous-time case and 

= /i(gi,g2,---,g7v), i = l,2,...,N 



(2.3a) 



(2.3b) 
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in the discrete-time case. Here, the tilde denotes the forward shift (m — t- m + 1) 
in the discrete-time coordinate m G Z. Note that the functions fi on the 
right-hand side of (12.31) are not ultralocal in their arguments qi, q2, ■ . . , Qn] 
they can be written in terms of {qi}, their partial derivatives or space/time 
shifts. In addition, in the (2 + l)-dimensional case, fi are, in general, non- 
local functions of {g,} involving a mixed action of d^, d~^, and dy, such as 

Let us start to construct the Miura maps from the Lax representation. 
Using the spatial part of the Lax representation fl2.2al) that is linear in ^ 
we obtain a set of coupled Riccati-type equations for \E'j\E'~\ i.e., 

I I 

k=l k=l 
= UiJ + Uu^^^J' -^^^J'U,,+ J2 Uik^k^^f 

- ^,^7! ^ f/,A*7\ l<ty^j<l (2.4a) 
in the continuous-space case (cf. [33]) and 

= Ui, + f/..^i,„^,-' + U^k'^k,n'^J:n, 1 < ^ 7^ J < ^ (2.4b) 

in the discrete-space case. The latter set of relations is derived from the triv- 
ial identity '^i^n+i'^J,n+i'^j,n+i'^J,l = '^i,n+i'^J^n- Here, (12.41) can be viewed as 
a linear algebraic system for the matrix elements Uij. Note that not all of the 
ratios between different components, ^i^J^ for i j, are independent quan- 
tities. We can suitably choose a set of Z — 1 distinct ^t^J^ (1 < 7^ j < 
so that all the other ^i^J^ can be expressed as a product of these quanti- 
ties and their inverse. In this section, we express this basis set of ^i^J^ as 
{pi,P2, ■ ■ ■ ,Pi-i}- Typically, we set pj := ^j+i^r^ pj := Thus, 
(12.41) can be reformulated as a system for {pi,p2, • • • , Pi-i}, their x-derivatives 
or spatial shifts, and {qi,q2, ■ ■ ■ ,qN}- Therefore, it is not ultralocal in the 
Pj, but it is ultralocal in the dynamical variables qj because the spatial Lax 
matrix U is assumed to be ultralocal in these variables. In this way, (12. 4p can 
provide a (typically) algebraic system of / — 1 equations for the unknowns 
qj. When I — 1 > N, one might consider that this system is overdetermined 
and has no solution. However, this is not the case as long as the original 
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Lax representation is consistent and provides a meaningful integrable system 
through the zero- curvature condition. In that case, the remaining Z — 1 — 
equations represent relations between some pj, for example, pj-^ is equal to a 
certain differential polynomial of pj^ . 

Similarly to the spatial part, the time part of the Lax representation 
(]2.2bp provides another set of coupled Riccati-type equations for ^fj\E'~^, i.e., 

for the derivation T2 in the continuous-time case and 

% {^^^-') = ( V,, + E V,u^,^.' [ V,, + E V,u^,^-' 

l<i7^j<l (2.5b) 

for the shift operator T2 in the discrete-time case. We can suitably take 
out / — 1 equations from (12.51) and reformulate them as the time evolution 
equations for {pi,P2, ■ ■ ■ ,Pi-i}- Note that these equations are not ultralocal 
in qj (cf. (11.3bp ) and depend on factors such as their partial derivatives or 
shifts, through the elements of the temporal Lax matrix V. 

Now, we have all the ingredients needed for constructing the (inverse) 
Miura maps: 

(1) evolution equations for the original variables {gi, q2, . . . , Qn} given 
by (O, 

(2) (typically algebraic) / — 1 "ultralocal" equations for {gi, ^2, • • • , Qn} ob- 
tained from (12.41) . wherein the coefficients involve {pi,p2, . . . ,pi-i}, 
their ^-derivatives or spatial shifts, 

(3) l — l evolution equations for the new variables {pi,p2, ■ ■ ■ ,pi-i} ob- 
tained from (12. 5p . which also involve {gi, g2, • • • , Qn}, their partial deriva- 
tives, shifts, etc. 

The key observation to be made (cf. (I1.3ap ) is that we can solve (2) to express 
min(/ — 1, A^) of the variables {qi, q2, . . . , qN} in terms of the remaining 
max(A^ — / -|- 1, 0) qj and min(Z — 1, A^) of the / — 1 variables {pi,P2, ■ ■ ■ ,Pi-i}- 
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These min(/ — 1, A^) expressions define a Miura map to the original system of 
evolution equations (1). The corresponding modified system is obtained 
by combining the max(A^ — / + 1, 0) evolution equations from (1) and the 
min(/ — l,iV) evolution equations from (3), wherein the unnecessary vari- 
ables must be eliminated using (2) in order to write the modified system in 
closed form. Note that (2) and (3) involve the spectral parameter (, so that 
the Miura map and the modified system also contain ( as an additional pa- 
rameter. The (nonstandard) Lax representation for the modified system can 
be obtained directly by substituting the defining expressions for the Miura 
map into the original Lax representation for (1); however, the free parameter 
C in the Miura map and the modified system must be distinguished from the 
spectral parameter ( in the original Lax representation. This is a brief sketch 
of our method; for clarity, we discuss the three cases / — 1 > A^, Z — 1 = A^, 
and / — 1 < A^ separately in more detail. 

• The case of / — 1 > A^ (Lax representation is "sparse" in the dependent 
variables): with a suitable renumbering of the pj, the Miura map is 
given by 

(pi,P2, ...,pn)^ iqi,q2, ■ ■ ■ ,g7v). 

Thus, the entire set of dependent variables is changed. In solving sys- 
tem (2) for a general value of (, we need to identify a basis set of A^ new 
variables pj that can be used to express the remaining I ~ 1 — N pj and 
the A^ original variables qj. When this elimination process turns out 
to be too complicated, we consider if fixing the parameter ( at some 
special value, say C = 0, can simplify system (2) to yield a parameter- 
less Miura map in concise form. In that case, the Miura map and the 
modified system should be written in terms of just A^ new variables pj, 
but these A^ variables are not required to form a basis set to express 
the other pj explicitly (cf. [SUES])- 

• The case of / — 1 = A^: the Miura map is given by 

{pi,P2,...,pi-i) ^ (gi,g2,...,g7v). 

Thus, the entire set of dependent variables is fully replaced without 
any need to choose a subset of the original/new dependent variables; 
conventionally, this is the case that the term "Miura map" refers to. We 
first solve system (2) to express {gi, q2, . . . , Qn} in terms of {pi,P2, ■ ■ ■ ,Pi~i}- 
Subsequently, substituting these expressions into (3), we obtain the 
modified system for the new variables pj. In the simplest case of 
/ = 2 and A^ = 1, that is, a 2 x 2 matrix Lax representation contain- 
ing only one dependent variable, this procedure is equivalent to Chen's 
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method [TSUin]. Note that the idea of Chen's method naturally comes 
from the original Miura map between the Gardner equation and the 
KdV equation |11[5] (also see Kruskal's article [57]). Chen's method 
can be suitably modified and applied to the discrete case [SHE^, the 
(2 + l)-dimensional case |60], and the multicomponent case |6T] as well. 

• The case of / — 1 < (Lax representation is "dense" in the dependent 
variables): with a suitable renumbering of the g^, the Miura map is 
given by 

(gi,g2, • • . ,gAr-i+i,pi,P2, • • • ,p«-i) ^ (gi,g2, • • • ,g7v)- 

Thus, only a subset of dependent variables is replaced. This appears 
to be the most interesting case and we concentrate on this case in 
the subsequent sections. System (2) of / — 1 equations is underdeter- 
mined, so that we can eliminate only / — 1 of the N original variables 
{qi, q2, . . . , Qn}, which are relabelled as {qN_i+2, ■ ■ ■ , Qn} above. That 
is, we mix — Z + 1 original variables qj with the I — 1 new variables 
Pj to form a closed system; actually, the way to eliminate I — 1 qj is, 
in general, nonunique. There exist maximally d^-^) different ways to 
eliminate I — 1 qj and mix the original and new variables. If the original 
system is totally asymmetric with respect to permutations of the vari- 
ables {qi, q2, . . . , qN}, then different eliminations may result in entirely 
different modified systems. Moreover, before applying the method, we 
can consider any invertible ultralocal change of the dependent variables 
in the original system; this change further leads to distinct modified 
systems. Therefore, we need to have an eye for beauty to sort out 
particularly interesting modified systems. 

In the first and second cases, / — 1 > A^, our method is conceptually similar 
to the method used in |62] (also see §3.1 of [63]). However, in contrast to our 
method, the method in [62|l63] uses components of the linear eigenfunction 
directly as the new dependent variables. 

Before ending this section, we remark that the definition of the new vari- 
ables {pi,P2, ■ ■ ■ ,Pi~i} intrinsically allows the following two freedoms: 

• the way of choosing p^ := \E'jj.\E'~^, which are algebraically independent 
and form a basis set to express all ^i^J^ (1 < 7^ J < 0' is nonunique; 

• for the Lax representation fl2.ip . an arbitrary gauge transformation 
ij) ^ (f) := gip can be considered, wherein g is a. constant invertible 
matrix to maintain the ultralocality of the spatial Lax matrix in the 
dependent variables. 
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Thus, these two freedoms represent a certain group of point transformations 
that act on the set of new dependent variables {pi,p2, . . . ,pi-i}. In the 
first case, / — 1 > A^, these freedoms may play a critical role in obtaining 
a simple modified system in closed form. In the second and third cases, 
/ — 1 < A^, they are not as essential as in the first case, because the entire set 
{pi,P2, ■ ■ ■ ,Pi-i} is used in the modified system. Nevertheless, a good choice 
of the new variables pj is crucial for obtaining a highly attractive modified 
system. 



3 (1 + l)-dimensional PDEs 

In this section, we apply the method described in section 2 to integrable PDEs 
in 1 + 1 dimensions. Four illustrative examples, namely, the Zakharov-Ito 
system, the Jaulent-Miodek system, the three-wave interaction system, and 
the Yajima-Oikawa system, are considered. 

3.1 The Zakharov— Ito system 

A two-component generalization of the KdV equation, 

Ut = Uxxx + 'iuUx + SwWx, (3.1a) 
wt = (uw)^, (3.1b) 

was proposed by Ito [35] in 1982. Actually, through the simple change of 
dependent variables q := —u/2 and r := —3w^/16, f l3.ip can be rewritten as 

qt = qxxx - Gqqx + 4:r^., (3.2a) 
n = -^xT - 2qr^, (3.2b) 

which was obtained earlier by Zakharov [13J. Thus, we may call fl3.ip or fl3.2p 
the Zakharov-Ito system. The Lax representation for (13. 2p is given by the 
pair of linear equations for a scalar function ^p [T3l|6l] , 

^xx={C + q + C^r) V^, (3.3a) 
^Pt = (4C - 2g)^, + q,^p. (3.3b) 

Indeed, the compatibility condition ip^xt = iptxx for (13. 3 p provides (13. 2p . We 
introduce a new variable v as ipx/'4' =• 'V- Then, the Lax representation (13. 3p 
implies the pair of relations 

v^ + v'^ = C + q + C^r, (3.4a) 
Vt = m-Hv + qx]x. (3.4b) 
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Relation fl3.4ap can define two different Miura maps, (g, v) (-)■ (g, r) and 
(f , r) f— )■ (g, r), both of which retain polynomiahty of the equations of motion. 

First, we consider the Miura map (g, v) i— i- (g, r) with r = —('^ + ({—q + Vx + f ^) 
(cf. |65] for the C = 1 case). Substituting this expression for r into fl3.2ap and 
combining it with ( 13.4b|) . we obtain the modified system for (g, f), 

qt = qxxx - Gqqx - + 4:({vxx + 2ffx), (3.5a) 
Vt = 4:(vx - 2{qv)x + qxx- (3.5b) 

Note that the hnear terms with the first-order x-differentiation can be ehm- 
inated using a Gahlean transformation. Indeed, by setting q ='■ q' — 2^, 
dt — 8(dx =: df, and dx =: dx', (13. 5p is simphfied to 

qt = qxxx - Gqqx + 4:C{vxx + '2vVx), (3.6a) 
Vt = qxx - 2(gw)x, (3.6b) 

where the prime is omitted for brevity. 

Second, we consider the Miura map {v, r) i— (g, r) with g = —( — (^^r + Vx + v"^. 
Substituting this expression for g into (13.4bp and (]3.2bl) . we obtain the mod- 
ified system for (f,r), 

Vt = QCVx + Vxxx - Qv'^Vx + 2C^{vr)x - C^Txx, 

n = 2(rx + QC'^rVx - 4:{vx + v'^)xr - 2{vx + v'^Yx- 
By setting r =: Cp^^ we can rewrite this system in the form 

Vt = 6Cvx + Vxxx - Gv^Vx + 2{vp^)x - {p^)xx, (3.7a) 
Pt = 2(px + Gp^Px - 2{vxP + v'^p)x- (3.7b) 

This two-component mKdV system, as well as the Miura map to the Zakharov- 
Ito system, is presented in section 4.2.17 of [66], wherein one can also find 
information on the C ~^ limit. A rather similar system is studied in [65] . 

3.2 The Jaulent— Miodek system 

Another two-component generalization of the KdV hierarchy was proposed 
by Jaulent and Miodek p6]. Its first nontrivial fiow is 

qt = rxxx + 2qxr + Aqrx, (3.8a) 
n = Aqx + Qrrx, (3.8b) 
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while the next flow allows the reduction r = to the KdV equation [30] ■ The 
Lax representation for fl3.8l) is given by the pair of linear equations, 



^xx + (g + Cr)V' = C'^, (3.9a) 
i^t = (4C + 2r)ij,, - r^i;. (3.9b) 

We introduce a new variable u as ipx/'^ ='■ u. Thus, the Lax representation 
(13. 9p implies the pair of relations 

M^ + M^ + g + Cr = C^ (3.10a) 
Ut=[{AC + 2r)u-r^]^. (3.10b) 

Relation f l3.10ap can define two different Miura maps to the Jaulent-Miodek 
system (13. 8p . i.e., (g, u) h- )■ (g, r) and (m, r) i— )■ (g, r); both of them retain poly- 
nomiality of the equations of motion. 

First, we consider the Miura map (g, u) t— )■ (g, r) with r = C, — C^^{q + + u^) 
Substituting this expression for r into (13. Sal) and (]3.10bl) and rescaling the 
time variable as dt =: — we obtain the modified system for {q,u), 

Qt Qx ~l~ Qxxx ~l~ QQQx ~I~ Uxxxx ~I~ (^ ^xxx 

+ 2g^M^ + 4gMa.a. + 2g^M^ + Squu^, (3.11a) 
U£ = -QC'^Ux - Uxxx + Qu'^Ux - qxx + 2{qu)x. (3.11b) 

Interestingly, (13. lip resembles a system of the coupled KdV-mKdV type. 

Second, we consider the Miura map (m, r) i— )■ (g, r) with q = ("^ — (r — — u"^ 
(cf. for ^ = and ^68j for general (). Substituting this expression for g 
into fl3.8bp and combining it with fl3.10bp . we obtain the modified system for 
(M,r), 

= 4Cmx - + 2(Mr)^, (3.12a) 
Vt = —4:(rx — -iuxx — 8uUx + Grr^. (3.12b) 

Note that the parameter ( in (I3.12p is nonessential and can be set equal to 
zero by the Galilean transformation r =: r' + 2(, dt — 8(8^ ='■ dt', ='■ d^'- 
Thus, (I3.12P is equivalent to the first nontrivial flow of the modified Jaulent- 
Miodek hierarchy studied in f6?] . 
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3.3 The three-wave interaction system 



We consider the three-wave interaction system [37j in the nonreduced form: 





- CiUi^x + (C2 


- C3)U3V2 


= 0, 


(3.13a) 




- C2U2,x + (C3 


- Ci)usVi 


= 0, 


(3.13b) 


U3,t 


- CsUs^x + (C2 


- Ci)uiU2 


= 0, 


(3.13c) 


Vl,t 


- CiVi^x + (C3 


- C2)V3U2 


= 0, 


(3.13d) 


V2,t 


- C2V2,x + (Cl 


- C3)V3UI 


= 0, 


(3.13e) 


V3,t 


- C3f 3,x + (Cl 


- C2)ViV2 


= 0. 


(3.13f) 



Here, the constants Cj {j = 1, 2, 3) are parametrized in terms of the constants 
aj and (3j as 

(3i-(32 

Cl = , C2 = , C3 = . (3.14) 

«! — 0^2 012 — 03 03 — ai 

Thus, they can be considered the "slopes" of the three sides of a triangle. 
Specific integrable systems describing triad wave interactions are derived 
by imposing a complex conjugacy reduction on f l3.13p to halve the number 
of dependent variables [121 169]. The Lax representation for the three- wave 
interaction system (13.131) with fl3.14p is given by the pair of linear PDEs |69] . 



^1 






Ml 


U3 




^1 


^2 






ia2C 


U2 




^2 


^3 _ 


X 


V3 


V2 






. ^3 


^1 " 






ClUi 


C3U3 




" ^1 


^2 




cm 


i/32C 


C2U2 




^2 


^3 . 


t 


_ C3V3 


C2V2 


i/33C . 




. ^3 



With ^2/^1 =: wi and ^3/^1 =: w^, (I3.15ap and fl3.15bp imply 

wi,x = vi + i(a2 - ai)Cwi + U2W3 - {uiWi + U3W2.)wi, 
W3,x = V3 + i(a3 - ai)(w3 + V2W1 - {uiWi + 1x3^3)^3, 



and 



wi^t = civi + i(/32 - f3i)(wi + C2U2W3 - (ciMiWi + 03^3^3)^1, 

W3,t = C3M3 + i(/33 - /3l)Cw3 + C2V2W1 - (CiUiWi + 03^3^3)^3, 



(3.15a) 



(3.15b) 



(3.16a) 
(3.16b) 



(3.17a) 
(3.17b) 
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respectively. Thus, relations f l3.16p can define the Miura map {ui, U2, M3, Wi, V2, w^) 

M2, U3, f 1, f2, "^3), which replaces two of the six dependent variables. 
Using fl3.16p . we can eliminate vi and V3 in fl3.17p to obtain 

wi^t - CiWi,^ + (ci - C2)u2Ws - (ci - Cs)u3WsWi = 0, (3.18a) 

W3,t - C3ty3,x + (C3 - C2)V2WI - (C3 - Ci)UiWiW3 = 0. (3.18b) 

In the same manner, (13.13bp and fl3.13ep can be rewritten as 

U2,t - C2U2,x + (Cs - Ci)u3 [wi^^ - 1(^2 - aOC^l - U2W3 + {uiWi + ^3^3)^1] = 0, 

(3.19a) 

V2,t - C2V2,x + (Cl - C3) [Ws^x - i(a3 - ai)C^^3 - V2W1 + {uiWi + U3W3,)W3] Ui = 0. 

(3.19b) 

The six equations, fl3.13ap . (13.1 3cl) . (13.181) . and (13.191) . comprise the modified 
system for {ui,U2, u^, Wi,V2, W3) in closed form. However, the system appears 
asymmetric and does not allow a complex conjugacy reduction directly. To 
rewrite the modified system in a more symmetric form, we apply a point 
transformation to some of the variables that were unchanged in the Miura 
map. 

Comparing (I3.13ap and (I3.13cl) with (I3.18ap and fl3.18bp . respectively, we 
can naturally move from the old variable f 2 to a new variable W2 as 

(C2 - 03)^2 =: (C2 - Ci)w2 + (Ci - C3)uiW3. (3.20) 

Here, we assume that Ci, C2, and C3 are pairwise distinct. Noting the identity 
(«! — a2)(c2 — Cl) = («! — 0:3) (c2 — C3), we can express the modified system 
for {ui, U2, U3, wi, W2, W3) in the desired form, 

- CiUi^x + (C2 - Ci)u3W2 + (ci - C3)uiU3W3 = 0, (3.21a) 

U2,t - C2U2,x + (C3 - Ci)u3Wi^x + Ki^l " a2)(c3 - Ci)u3'Wi 

- (C3 - Ci)u2U3W3 + (C3 - Ci){UiWi + U3W3)u3Wi = 0, (3.21b) 

U3,t - CsUs^x + (C2 - Ci)uiU2 = 0, (3.21c) 
Wi,t - CiWi^x + (Cl - C2)W3U2 + (C3 - Ci)WiW3U3 = 0, (3.21d) 

W2,t - C2W2,x - (Cl - C3)w3Ui^x + iC(«l " a2)(Cl - 03)^3^1 

- (Cl - 03)1^2^3^*3 + (Cl - C3){uiWi + U3W3)W3UI = 0, (3.21e) 

W3,t - C3W3^x + (Cl - C2)wiW2 = 0. (3.21f) 

The "refined" Miura map {ui,U2,U3,Wi, 102,103) (mi, M2, M3, Vi, ^2, "^^3) from 
(K2T^ to fl313D is given by combining flHTTHj) and (KW^ . Actually, it is 
possible to eliminate the quartic terms in fl3.21bp and (I3.21ep by further 
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applying a nonlocal transformation. Indeed, the simplest conservation law 
for dSSI]), 

{uiWi + U3W3)t = (ciMiWi + CsUsWs)^, 

motivates us to introduce the new set of variables (gi, g2, Q's, ''"i, ^"2, '"a) by 

(3.22) 

Here, S is an arbitrary constant. Note that exp [f^{uiWi + U3Ws)dx'~\ can 
be identified with exp(— iai^x — i(3i(t) in the Lax representation fl3.15l) . 
up to a multiplicative constant. Under suitable boundary conditions, the 
modified system fl3.2ip is transformed to the form 

qi,t - cigi.x + (c2 - Ci)g3r2 = 0, (3.23a) 
q2,t - C2q2,x + (c3 - Ci)g3ri,^. + iC(ai - 0^2) (cs - c^ggri 

+ 6{ci- 02)^1^2^1 + [(ci - C3) + 5(c3 - C2)] g2g3^3 = 0, (3.23b) 

gs,* - csgs,!^ + (c2 - ci)gig2 + (1 + 5)(ci - 03)^1 gsri = 0, (3.23c) 

ri,t - ciri,^ + (ci - C2)r3g2 = 0, (3.23d) 
r2,t - C2r2,x - (ci - c^y^qi^^ + i({ai - 02) (ci - Csy^qi 

+ 5(c2 - ci)rir2gi + [(03 - ci) + 5(c2 - C3)] r2r3g3 = 0, (3.23e) 

r3,t - c^r^^^ + (ci - C2)rir2 + (1 + 5)(c3 - Ci)rir3gi = 0. (3.23f) 

Some special cases, such as 5 = —1 or and C = 0, appear to be particularly 
interesting. 

Note that fl3.20p is not the only point transformation that can convert the 
modified system for {ui, U2, us, wi, V2, W3) to a symmetric form that allows a 
complex conjugacy reduction. Let us consider the simplest case of C = and 
change the variables U2 and f 2 to U2 and W2 as 

^ ai — ^ 

U2 - U3W1 =: U2, V2 - U1W3 =: W2. (3.24) 

ai — a2 
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Thus, the modified system for {ui,U2,Us,Wi,W2,W3) can be written as 

(3.25a) 



Ul,t - Citii.x + (C2 - Ci)u3W2 + (C2 - C3)UiUsW3 = 0, 
U2,t - C2U2,x + (Cl - C2)uiU2Wi + (C2 - C3)u2U3W3 

+ (C3 - C2) [(u3Wl)x + {uiWi - U3W3)u3Wi] = 0, 

U3,t - C3?^3,x + (C2 - Ci)UiU2 + (C2 - Ci)UiU3Wi = 0, 
- CiWi,^ + (Ci - 02)^3^2 + (C3 - C2)wiW3M3 = 0, 

W2,t - C2W2,x + (C2 - Ci)WiW2Ui + (C3 - C2)W2W3U3 
+ (C3 - C2) [(U^3^^l)x - (^ilW^i - U3W3)w3Ui] = 0, 

W3,t - C3W3,^ + (Ci - C2)WiW2 + (Ci - 02)^1^3^1 = 0. 



(3.25b) 
(3.25c) 
(3.25d) 

(3.25e) 
(3.25f) 



The associated spectral problem in the canonical form can be obtained by 
applying a gauge transformation to (13.1 5al) with (I3.16P and (13.241) . It is given 
by 



" $1 " 




iaiC + "^^i"^! + U3W3 


(^2 - Oii)Ui 








ia2C ~ UiWi 


. "^3 . 


X 


KW3 


W2 



02— ai 



U2 



i«3C - U3W3 





' $1 " 




$2 







(3.26) 



Thus, (13.241) is a natural point transformation from the point of view of a 
Lax representation. In the same way as for (13.211) . we can also consider a 
nonlocal transformation like (I3.22p . Indeed, using 



-7 (uiwi+u^w^)dx' 

1 

7 (Miuii+M3Ui3)da;' 



U2e 



-S {u-iW\+u^W'i)dx' 



W2e 



(7+'5) (uiwi+U3W3)dx' 



as the new set of variables and choosing the parameters 7 and 6 appropriately, 
we obtain a simplified version of the modified system (I3.25p . In particular, 
it is possible to eliminate the cubic terms in (I3.25ap . (I3.25cp . ( I3.25dp . and 

(HSMj). 



3.4 The Yajima— Oikawa system 

We consider the Yajima-Oikawa system |38] written in the general nonre- 
duced form, 



iQt + Qxx -PQ = 0, 
iRt — Rxx + RP = O, 
iPt + 2{QR), = 0. 



(3.27a) 
(3.27b) 
(3.27c) 
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The Lax representation for (I3.27p is given by the set of hnear equations 

(cf. HO]), 

= C^ + P^ + Q$, (3.28a) 

= m, (3.28b) 

i^t + Q$ = O, (3.28c) 

i<l>t + - R^^ - = O. (3.28d) 

Because the hnear eigenfunction comprises two components \1/ and $, there 
exist two distinct ways to define a Miura map to f l3.27p . 

First, we set =: R and '^x^~^ ='■ P- Thus, the spatial part of the 

Lax representation, ( ]3.28bp and fl3.28ap . imphes that 

R = Rx + RP, P = -CI + P^ + p2 - QR. 

These relations define the Miura map {Q, R, P) i— )■ {Q,R,P). Using fl3.28ap 
and f l3.28bp . the time part of the Lax representation, f l3.28cp and (]3.28dp . 
provides the evolution equations for R and P. Combining them with f l3.27ap . 
we arrive at the modified system for {Q, R, P), 

iQt + CQ + Qxx - PxQ - P^Q + QRQ = O, (3.29a) 
iRt - C-R - Rxx - RPx + RP"^ - RQR = O, (3.29b) 
iPt + {QR)x + QRP - PQR = O. (3.29c) 

In the case of scalar variables, this system with C = was studied in [7HI72]. 
and the Miura map to the Yajima-Oikawa system f l3.27p is also known |73] . 
Note that the parameter ( is nonessential if we consider the modified system 
fl3.29p separately. 

Moreover, we can transform fl3.29p to a simpler form. Indeed, the pair of 
relations = -P^ and i'^t = —QR"^ motivates us to introduce the new set 
of variables as 

q ■= ^"iQ, r := = $, p := ^"^P^ = ^"^^a.. 

Then, it is easy to show that they satisfy the closed system, 

m + + Qxx + '^PQx = O, 
in -Cr- Txx + 2r^.p = O, 
+ {<lf)x + Vl''^ ~ I'^V = o, 

which could be called the derivative Yajima-Oikawa system. Again, the 
scalar case with C = was studied in [7HI72]. 
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Q. Because (I3.28bp implies the relation 



Second, we set ^$"i(= R'^) 
$^$-1 = RQ^ f lX^ gives 



Q = Qxx + 2Q^RQ + Q(RQ) + QRQRQ - CQ - PQ- 



(3.30) 



This relation defines the Miura map {Q,R,P) i— )■ {Q,R,P), which changes 
only one of the three variables. Using also the time part of the Lax repre- 
sentation, fl3.28cp and (I3.28dp . we obtain the evolution equation for Q, 



iQt + Q.. + 2(QRj Q-PQ = 0. (3.31) 
Substituting fl3.30p into f l3.27cp . we obtain the evolution equation for P, 



iPt-2C[QR] +2[Q,,R 



+ 



4:Q^RQR + 2Q (RQ) R-2PQR + 2(QR 



O. 



(3.32) 



The three equations (13.311) . (I3.27bp . and (13.321) comprise the modified system 
for {Q,R,P). However, this modified system has no physical significance in 
its present form, because it looks asymmetric with respect to Q and R. That 
is, we cannot directly impose a conjugate relation between the two variables. 
To restore a desired symmetry, we only have to change the variable P as 
P — 2QxR + (I =: P, where the tilde does not denote the forward shift in 
the discrete-time case. Thus, the modified system takes the form 



iQt + CQ + Qxx + 2QR^Q -PQ = 0, 
iRt -CR- Rxx + 2RQ^R + RP = 0, 
iPt + 2(Q,R^ -2(qR,] QR + 2QR[Q,R 



2PQR, - 2Q^RP + 2 



QR 



O. 



(3.33a) 
(3.33b) 

2QxRQxR + 2QRxQRx 
(3.33c) 



The Miura map to the original Yajima-Oikawa system (I3.27p is given by the 
relations 

Q = Qxx + Q [rq) ^ + QRQRQ -PQ, P = -QI + 2g,i? + P. 

Note that by setting P = P' + a{Q^R - QR^) + /3{QRf, we obtain a two- 
parameter deformation of the modified system (13.331) . 
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4 Differential- difference equations 



In this section, we apply our method to differential-difference equations in 
1 + 1 dimensions, such as the Toda lattice, the Belov-Chaltikian lattice, 
the relativistic Toda lattice, and the Ablowitz-Ladik lattice, wherein the 
spatial variable n takes discrete values. We can also apply our method in 
the fully discrete case, wherein both the space and time variables are dis- 
cretized. However, the results appear to be too complicated compared with 
the continuous-time case, so we do not discuss the fully discrete case in this 
paper. 

4.1 The Toda lattice in Flaschka— Manakov coordinates 



We consider the Toda lattice written in Flaschka-Manakov coordinates [391 - 
E]: 

Un,t = Un{Vn- Vn-l), (4.1a) 
Vn,t = Un+1 - Un- (4.1b) 

The parametrization m„ = e^"~^"-i, Vn = Xn,t enables the system (14.11) to be 
rewritten as the Newtonian equations of motion for the Toda lattice. 

The Lax representation for the Toda lattice (14.11) [TO] can be written in the 
scalar (or "big matrix") form, 

^^^+1 + u„^„_i = (C + i;n)*n, (4.2a) 

'^n,t = (4.2b) 

Indeed, the commutativity of the spatial shift and the time derivative for \E'„ 
results in (14. ip . We introduce a new variable g„ as '^n+i/'^n ='■ —Qn- Then, 
the Lax representation fl4.2p implies the pair of relations 

Un = + Vn)qn-i - QuQu-i, (4.3a) 

qn,t = -Un+1 - Qn [(C + ^n) + Qn] ■ (4.3b) 

The first relation fl4.3ap defines the Miura map (g„,f„) ^ {un,Vn)- Using 
(I4.3ap . we can eliminate and u„ in fl4.3bp and fl4.1bl) to obtain a closed 
differential-difference system for (g„,t>„), 

qn,t = qniqn+l+Vn+l- Qn-Vn), (4.4a) 
(g„+l + Vn+l + C) + qn-l{qn + Vn + () ■ (4.4b) 
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Note that the parameter ( is no longer essential if we consider f l4.4l) separately 
as an isolated system. Indeed, it can be set equal to zero by shifting t>„ as 
Vn + C ='■ 'v'n- III addition, it is easy to identify (14 ■4p with the Volterra lattice 
if we use the new pair of variables {—qn,(ln + Vn + C)- 

The two-component system (14.11) allows a one-parameter generalization 
corresponding to the relativistic deformation of the Toda lattice (see Chapter 
6 of [in])- Thus, it would be interesting to obtain the one-parameter gener- 
alization of (14.41) along the same lines as above and discuss its relationship 
with the relativistic Volterra lattice |10j . 



4.2 The Belov-Chaltikian lattice 

The Belov-Chaltikian lattice |42j is described by the equations of motion 



Un,t = Un{Un+l - U^^l) - (w„+l - Wn), (4.5a) 
Wn,t = Wn{u„,+ 1 - Un-2)- (4.5b) 

The Lax representation for the Belov-Chaltikian lattice (14.51) is given by the 
pair of linear equations (cf. [74j). 



- M„^„ + Wn'^n-l = (4.6a) 
^n,t = + Un-l^n- (4.6b) 

If we introduce a new variable f„ as Vn '■= '^n+i/'^n, the Lax representation 
(14.61) implies the pair of relations 

U)n = {(Vn+lVn - Vn + Un)Vn-l, (4.7a) 
Vn,t = Vn{Vn+l + Un) - fn(fn + Mn-l)- (4.7b) 

The first relation (I4.7al) defines the Miura map {un^Vn) [un^Wn)- Using 
(I4.7ap . we can eliminate Wn+i and w„ in (I4.5ap : combining it with (]4.7bp . we 
arrive at a closed differential-difference system for («„,!>„), 

Un,t = Un{Un+l " Un-l) " ^n+l^^n + Mn^n-l 

+ Vn{Vn+l - Vn-1 - CVn+2Vn+l + C^^n+l^^n-l) , (4.8a) 

Vn,t = Vn{Un - ^n-l + - Vn). (4.8b) 

In the special case = 0, (I4.7ap coincides with the Miura map proposed 
in |75| (also see [76]) for the Bogoyavlensky lattice, although the reduction 
M„ = is not consistent with this particular flow. 



23 



4.3 The relativistic Toda lattice 

The relativistic Toda lattice introduced by Ruijsenaars 



g'2 qX^^i Xjl 
1 + (y(2e^n + l- 



•^n,t-^n—l,t 



1 + g^e 



n Xji — l 



can be naturally written in the two-component form 



e''"(l + /e^"+^-^" 



7-) — n^ci^ri-\-l Xn-\-pn 



^2^X77, Xn—X -\-pn—l 



The 2x2 Lax representation for fl4.1Up is given by [77] 



g^Q~Xn'-\-'pn Q 



(^-2 _|_ ^2gX„-x„_i+p„_i _(^-lgi.'»i 
^-1^2g-a:„_i+p„_i g 



^2,n 



(4.9) 



(4.10a) 
(4.10b) 



(4.11a) 
(4.11b) 



We introduce a new variable i/n in the exponential form e^" := ^l'i,n/^2,n 
so that a symmetric structure in the new pair of variables Xn and y„ can be 
uncovered. The Lax representation (14. lip enables p„ and yn,t to be expressed 
as 



^-1 _ 



n yn 



_|_ g2Q-X„+yn + l 



yn,t = C +ge 



(4.12a) 
(4.12b) 



Thus, fl4.12al) defines the Miura map {xn,yn) ^ ixn,Pn)- Using fl4.12ap . we 
can eliminate pn in fl4.10ap and Pn-i in (I4.12bl) to obtain a closed differential- 
difference system for {xn,yn), 



yn,t 



_j_ g^^yn+l—Xn 

{(-^ - e^"-2'")(l + ^2gj/„-y„_i^ 

^ _|_ g^^y-n Xyi — \ 



(4.13a) 
(4.13b) 



Actually, ( I4.13P gives an auto-Backlund transformation between the two so- 
lutions Xn and Un of the relativistic Toda lattice (14. 9p . 
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4.4 The Ablowitz-Ladik lattice 



In this subsection, we consider the Ablowitz-Ladik lattice [H], which appears 
to be the most instructive example in the discrete-space case. The equations 
of motion for the (nonreduced form of the) Ablowitz-Ladik lattice are 



Qn,t — O-Qn+l + bQn-1 + (a — b)Qn + aQn+lRnQn — bQnRnQn-1 
Rn,t — bRn+1 + aRn-l + {b — a)Rn + bRn+lQnRn — CtRnQnRn-1 " 

and the Lax representation is given by 



= 0, 
(4.14a) 

--0, 
(4.14b) 



^l,n+l 



Ch Qr. 

Rn tI; 



^2,n 



(^2 _ _ aQ^R,^ 



CaQn + jQ 



n-l 



(aRn-l + jRn (72 — 1)^2 — bRnQn-l 



(4.15a) 
(4.15b) 



Here, a and b are free parameters. Note that the Ablowitz-Ladik lattice 
fl4.14p is integrable for matrix-valued dependent variables [78] (also see [79|l80] 
and references therein); in the general case, Qn and i?„ are li x I2 and I2 x li 
matrices, respectively. We consider an (/i + I2) ^ h matrix- valued solution 
to the pair of linear equations (14.151) such that \E'i,„ is an li x li invertible 
matrix. Then, in terms of the I2 x li matrix Pn := '^2,n^ili-, (I4.15P can be 
rewritten as a pair of discrete and continuous matrix Riccati equations for 



Rn 



n+1 



c 



P-n ~l~ Pn+lQnPn 



(4.16a) 



ra,t 



CaRn-l + -Rn - - l)aP„, + {-:-l]bPn 



c 



- bRnQn-lPn + aPnQnRn-l - -^PnQn-lPn " (aPnQnPn- (4.16b) 

The first relation fl4.16ap defines the Miura map {Qn, Pn) ^ (Qn, Rn)- Using 
fl4.16ap . we can eliminate Rn and Rn-i in (14. Map and (]4.16bp to obtain a 
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closed differential-difference system for {Qn,Pn) as (see (3.16) in 

Qn,t - aQn+l + hQn-l + (o - h)Qn + aQn+l \CPn+l " \Pn \ Qn 



c 

^Qn I CPn+l ~ -jPfi 1 Qn-1 + O'Qn+lPn+lQnPnQn ~ bQnPn+lQnPnQn-1 = O, 



c 



i'n.t - bPn+1 + aP„_i + (6 - a)P„ - bPn+l { ]:Qn ' CQn-1 ] Pn 



(4.17a) 



c 



OiPn \ ^Qn CQn—l 1 P-n—l ~l" ^-Pn+lQn-PnQn— l-^n OiPnQ nPnQ n—lP-n—l O. 



c 



(4.17b) 



When a = — 6 = i, fl4.17p gives an integrable space discretization of the Gerdjikov- 
Ivanov system (11.41) . 

Note that the parameter C, in the modified system (14.171) is nonessential 
as in the continuous case. Indeed, for the elementary cases a = or 6 = 0, a 
transformation of the form 

Qn = e^e^'Q'n. Pn = C^""^' e'^' P^ 

where c is a suitably chosen parameter, and a rescaling of t can be used to 
set C as 1. For the more general case of a6 7^ 0, we consider the redefinition 
of the parameters C^a = a', = b', instead of rescaling t. 

5 (2 + l)-dimensional PDEs 

In this section, we illustrate how to apply our method in 2 + 1 dimensions. 
As instructive examples, two distinct generalizations of the NLS system to 
2 + 1 dimensions are considered. 

5.1 (2 + l)-dimensional NLS: Calogero— Degasperis sys- 
tem 

In their pioneering paper [IBj , Calogero and Degasperis proposed a large class 
of multidimensional PDEs that can be associated with the same spectral 
problem as the (1 + l)-dimensional NLS system. A representative example 
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of the class is a (2 + l)-dimensional generalization of the NLS system fll.ip . 



iQt + Qxy - fQ - Qg = O, 
iRt - R^y + gR + Rf = 0, 



fx 

9x 



iQR)y, 
{RQ)y. 



(5.1a) 
(5.1b) 
(5.1c) 
(5.1d) 



Here, we consider the general case where the dependent variables are matrix- 
valued (cf. [51]). Using fl5.1cp and (I5.1d|) . the auxiliary fields / and g can be 
formally written in the nonlocal forms / = d~^dy{QR) and g = d~^dy{RQ), 
which can be substituted back into fl5.1al) and (]5.1bl) . Note that by setting 
t = t„+i and y = tn [HI], (15.11) can be identified with the recursion relation 
for the (1 + l)-dimensional NLS hierarchy [53]. Conversely, any (1 + 1)- 
dimensional integrable system having a proper recursion operator can be 
generalized to 2 + 1 dimensions in this way. 

In the literature, the Calogero-Degasperis system (15.11) is often referred 
to as the Zakharov system [IB]. The Lax representation for this system is 
given by [ISlEniES] (cf. O) 



^1 

^2 

^1 

^2 



-iC/i Q 
R iCh 



2C 



^1 

^2 



+ 



^1 

-if iQy 
-iRy ig 



^1 

^2 



(5.2a) 
(5.2b) 



Here, the spectral parameter which is independent of x, has to satisfy the 
non-isospectral condition Q = '^CCy [22] . The Lax representation (15.21) is valid 
for the simpler isospectral case (t = Cy = 0- However, the non-isospectral 
nature of the parameter ( turns out to be explicit and essential in applying 
the inverse scattering method [IB]. 

In terms of the new variable P := \l/2^r^i Lax representation (15. 2p 
can be reformulated as a pair of matrix Riccati equations. 



Px 

Pf 



-- R + 2iCP - PQP, 
2CPy = -iRy + igP 



iPf - iPQyP 



(5.3a) 
(5.3b) 



Thus, using (I5.3ap . which actually defines the Miura map {Q,P) ^ {Q,R), 
we can eliminate R from (I5.1cp . fl5.1dl) . and fl5.3bp to obtain a closed system 
for [Q, P) with the auxiliary fields / and g. To express this system in a more 
symmetric fashion, we redefine the auxiliary fields as 



/ 



Cyh + l{QP)y + u, 



'^Cyh + l{PQ)y + V. 
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Then, we obtain a (2 + l)-dimensional generalization of the Gerdjikov-Ivanov 
system fll.4p in the form 



iQt + KyQ + Qxy - |w + liQP)y^ Q - Q + ^(Pg)j,| = O, 

iPt + KyP - p^y + !^v- ^(pg),} p + p |m - ^(gp),| = o, 



-2icgp + -(gp. - g.p) + {QPf 
-2icpg + i(Px.g - pg.) + {pqY 



0. 



(5.4a) 
(5.4b) 
(5.4c) 

(5.4d) 
(5.4e) 



It is now clear how to impose suitable reductions, such as the complex conju- 
gacy reduction, on the dependent variables. In the case of = Cy = 0, fl5.4ap 
and (15.4bp imply the conservation law 



-(QyP-QPy) 



d 
dy 



-{QxP-QPx)-{QPf 



0. 



Using the same technique as in the (1 + l)-dimensional case (see (11. 6p ). we 
can construct an infinite set of conservation laws for the Calogero-Degasperis 
system (15.11) (cf. p3j). Thus, the conservation laws for (15. 4p in the general 
case of nonconstant C, can be obtained with the aid of (I5.3ap . 



5.2 (2 + l)-dimensional NLS: Davey— Stewartson system 



The other (2 + l)-dimensional generalization of the NLS system (II. ip to be 
considered is (the integrable case of) the Davey-Stewartson system [18], also 
referred to as the Benney-Roskes system [81]; it can be written in the form 



iQt + aQxx + bQyy 

iRt — ClRxX ~ bRyy 

= 2b{QR)y, 

Qy = 2a(Pg)^. 



fQ-Qg = 0, 
gR + Rf = 0, 



(5.5a) 
(5.5b) 
(5.5c) 
(5.5d) 



Here, a and b are free parameters, and / and g are auxiliary fields. The sim- 
plest case of 6 = (or a = 0) provides a (2 + l)-dimensional generalization of 
the Yajima-Oikawa system (I3.27p [I31IB5]- Note that the Davey-Stewartson 
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system fl5.5p is integrable for matrix- valued dependent variables [SUHSHSZ]- 
Its Lax representation is given by P ^ ITO t lSS t lS^ 







' 


Q ' 




' ^1 " 


. ^.^2 . 




R 







^2 



" ^1 " 




^2 


t 



-hdl 



f 



-hRdy + hRy 



aQdx — aQx 
adl-g 



^1 

^2 



(5.6a) 
(5.6b) 



In contrast with fl5.2al) . the spatial part of the Lax representation, fl5.6al) . is 
a two-dimensional problem without a spectral parameter. 

Let us introduce a new variable P := ^!2^i^. Owing to fl5.6al) . we can 
introduce the auxiliary field w as 



(5.7) 



Note that {-^l 



QP. The compatibility condition of these linear 



PDEs for ^1 implies the relation 

Wx - {QP)y + [W, QP] = O. 

From fl5.6ap and fl5.7p . we obtain 

R = P^ + Pw, 



(5.8) 
(5.9) 



which defines the Miura map {Q,P) ^ {Q,R). Thus, equations fl5.5cp and 
(15.5dp for the auxiliary fields / and g can be rewritten as 

fx = 2biQPy + QPw)y, (5.10) 
gy = 2a{PyQ + PwQ)^. (5.11) 

Noting the identity 

^2,.x*r' = (^2^r')..-2(^2^r')x^i(^r')x+2^2(*r')x.^i(^r')x-^2(^r') 

and using (15.6bp . we can compute the time derivative of P{= '^2'^i^) with 
the aid of ([52D and ([52D as 

iPt = a^2,xx^r^ - ^*2^r^ - ^^*i,2/*r^ + bRy + 6^2^'r^^i,ra^r^ 
- ^2^rV - a^2^r^<5^2,x^r^ + a^2^r^Qx^2^r^ 

= aPxx + 2aPxQP + 2aPQPQP - aPQPQP + aP{QP)x 
-gP- bPyW - bPw^ + bPyy + b{Pw)y + bPwy + bPw^ 
-Pf- aPQPx - aPQPQP + aPQ^P 

= aPxx + bPyy -Pf-gP + 2bPwy + 2a{PQ)xP. (5.12) 
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Thus, fl33al) . (Eg), fICTJ]) . and f lETT]) together with (ESD comprise the mod- 
ified system. To rewrite it in a more symmetric fashion, we redefine the 
auxihary fields as 



Then, we obtain a (2 + l)-dimensional generahzation of the Gerdjikov-Ivanov 
system (11. 4p in the form 



iQt + aQxx + hQyy + Q[F- a{PQ),] + {H - hwy) Q = 0, (5.13a) 
iPt - aP,, - bPyy - [F + a(Pg) J P-P{H + bwy) = O, (5.13b) 




For the scalar (and thus commutative) case, this system and the Miura map 
to the Davey-Stewartson system were derived in [90] using a different ap- 
proach. 

6 Concluding remarks 

In this paper, we have developed an effective method of identifying new in- 
tegrable systems that can be mapped to a given integrable system by Miura 
transformations. The method is applicable in a systematic manner as long as 
the spatial part of the Lax representation for the given system is ultralocal in 
the dependent variables. In fact, most of the known integrable systems sat- 
isfy this requirement, and the wide applicability of the method is illustrated 
using numerous examples. The cornerstone of our method is to overcome a 
stereotype of perceiving the original dependent variables and an eigenfunc- 
tion of the associated linear problem as entirely different objects that cannot 
be swapped. That is, we combine a subset of the original dependent vari- 
ables and the components of an eigenfunction of the linear problem to form 
a new set of dependent variables that satisfies a closed system; this process 
elucidates the true nature of the Lax representation and explains its origin. 
In the appendices, a variant of the method is used to derive and characterize 
derivative NLS systems of the Chen-Lee-Liu type. 

To illustrate the method, we mainly considered the first nontrivial fiows of 
integrable hierarchies as examples. Note, however, that the method can also 
be applied to the higher fiows in each integrable hierarchy, as well as to the 
negative fiows. Indeed, a Miura map between two integrable systems actually 



/ =: bwy -H, g=: a{PQ), - F. 




Fy + a {PyQ - PQy + 2 P w Q) ^ = O, 

+ b {QPy - QyP + QPw + wQP)y = O 
- {QP)y + [w, QP] = O. 



(5.13c) 
(5.13d) 
(5.13e) 
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gives a link between the original hierarchy and the modified hierarchy rather 
than between a particular flow and its modification. Thus, an integrable 
system is not an isolated object either "longitudinally" or "transversely"; it 
not only arises as a member in an infinite hierarchy of commuting fiows, but 
also appears with its partners related by a chain of Miura maps. That is, 
a given integrable system is either an original system having one or more 
modifications or a modification of a more basic system or both. 

Originally invented to derive new integrable systems from known ones, 
our method also provides the most direct route to solving the derived inte- 
grable systems. By construction, the set of dependent variables satisfying a 
derived modified system is written explicitly in terms of the original variables 
and an eigenfunction of the linear problem for the original system. More- 
over, when the inverse scattering method is applied to the original nonlinear 
system, it intrinsically provides formulas for determining the eigenfunctions 
of the associated linear problem. Thus, if the original system can be solved 
by the inverse scattering method, then it is an easy and straightforward task 
to obtain a solution formula for the modified system using simple operations 
such as division and integration by parts. Incidentally, this kind of operation 
plays a critical role in integrating the sine-Gordon equation and its discretiza- 
tions explicitly. Note, however, that the general solution of an / x / matrix 
or /-th-order scalar linear problem is expressed as a linear combination of / 
fundamental solutions with I arbitrary coefficients. Therefore, we need to 
impose suitable boundary conditions on the modified system to remove this 
arbitrariness and obtain its solution formula uniquely. 

For a (1 -|- l)-dimensional system, the existence of an arbitrary param- 
eter, called the spectral parameter, in the Lax representation is a manifes- 
tation of its integrability. Consequently, the Miura map and the modified 
system obtained by our method naturally contain the spectral parameter 
as a free parameter. In general, the spectral parameter has its origin in 
a group of point transformations that leave the original unmodified hierar- 
chy invariant. For example, the spectral problem f ll.2ap is invariant under 
the one-parameter group of transformations: := \E'ie^'^'^', \E'2 := \E'2e~^'^^, 
Q' ■= Qe^i'^^, R' := Rq-'^'^'' , and C' := C - k. Thus, the associated NLS hi- 
erarchy is also invariant under the transformation of dependent variables, 
Q' = Qe^""^^ and R' = Re'"^^^^. Note that this transformation induces a lin- 
ear change in the infinite set of time variables for the NLS hierarchy, and 
thus, each fiow is not invariant. The invariance of individual NLS fiows can 
be restored by applying the inverse of this linear change in the time vari- 
ables. In the simplest case of the first nontrivial fiow, the NLS system fll.ll) . 
this amounts to its Galilean invariance [91j. This implies that the parameter 
( in the corresponding modified system flL4p can be fixed at any value by 



31 



a point transformation, cf. appendix \M In general, the spectral parameter 
appearing in the Miura map is nonessential in the sense that it can be fixed 
at any generic value using a point transformation. However, such a point 
transformation usually does not leave the modified integrable hierarchy in- 
variant, so the spectral parameter can be viewed as a deformation parameter 
of the modified hierarchy. 

To avoid needless confusion, we did not proactively address the issues of 
certain types of nonstandard Lax representations; in fact, there exist cases 
wherein the components of the linear eigenfunction in the matrix Lax repre- 
sentation are not entirely independent. Thus, we have to identify an appro- 
priate subset of the components to define the relevant Miura map correctly; 
this point was briefiy and partially touched upon in section [2J There are two 
main types of such nonstandard Lax representations. In the first type, linear 
ordinary differential/difference relations among components with respect to 
the spatial variable, such as \l/3 = hold true. Typical examples of this 

type are high-order scalar Lax pairs reformulated in the matrix form. In the 
second type, we have nonlinear algebraic relations among components, such 
as (^2)^ = ^1^3; so that not all components are algebraically independent. 
For example, the scalar case of the NLS system f lLlI) allows the standard 
Lax representation in the 2x2 matrix form (see (11.21) ): however, it also im- 
plies a 3 X 3 nonstandard Lax representation for the "squared" eigenfunction 
((^1)2,^1^2,(^2)') [inland a 4 X 4 nonstandard Lax representation for the 
"cubed" eigenfunction with components (\l/i)^~-'(\l/2)-', j = 0,1,2,3 and so 
forth. For a given Lax representation, there exists no algorithmic way to find 
all nontrivial relations among components of the eigenfunction. However, it 
is often possible to identify such relations by noting the internal symmetry 
of the Lax pair. For example, if the Lax matrices are antisymmetric up to 
a certain similarity transformation, with respect to the secondary diagonal, 
then a quadratic relation among the components results (cf. |61]l92]). 

As a final remark, we note that our method can also be applied to classi- 
cal many-body problems such as the Calogero-Moser models. In that case, 
the method pinpoints a nontrivial transformation from the modified system 
to a given original system. The transformation is no longer a Miura map in 
the usual sense of the word; rather, it is a coordinate transformation acting 
on the finite set of dynamical variables {q,p}- The existence of such a trans- 
formation is highly nontrivial; it maps the modified system to the original 
system in an easy-to-follow manner, but its inverse cannot be found without 
using a systematic approach. 
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A Continuous Chen— Lee— Liu system 

In this appendix, we show that the Lax representation for the NLS system 
can provide the solutions of a derivative NLS system, called the Chen-Lee- 
Liu system This result can be obtained by exploiting the transformation 
of dependent variables between the Gerdjikov-Ivanov system (11. 4p and the 
Chen-Lee-Liu system. 

Let Q and R satisfy the nonreduced NLS system (11.11) . Let and ^2 be 
the first and second components of an eigenfunction of the associated linear 
problem (II. 2p . Then, we have the following proposition. 

Proposition A.l. The new pair of variables, 

{q := [^,e<^^^'^'Y 
\ r := ^ae^^^+'^'^'S 

satisfies a closed system, i.e., the derivative NLS system 

ig* + Qxx + 4iCgrg + 2qrq^ = O, 
ift — r^x — 4:i(rqr + 2rxqr = O. 



(A.l) 



(A.2) 



Remark. The exponential factor e^^^"*"^^^ * is introduced in (lA.ip simply 
to remove nonessential linear terms from the resulting system, (lA.2p . The 
standard form of the Chen-Lee-Liu system corresponds to the case of C = 
(see [20H2I] for the matrix case). Note that the parameter ( is nonessential 
and flA.2p can be reduced to this case by a point transformation. Indeed, if 



we change the variables in (II. 2p as 

:= ^^e'^^+2iC'*, := ^ae'^'^^-^"^'*, Q' := Qe^'<''+^'<'\ R' := i^e-^'^^-^'^'*, 
df := dt - 4(9^, d^' := 

and omit the prime, we obtain the same Lax representation (II. 2p with C = 0. 
This is a manifestation of the Galilean invariance of the NLS system (II. ip |91] . 
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Proof. For brevity, we use the quantities $j := vl/^e'''^'^^'^^* (j = 1, 2) instead 
of Using the Lax representation fll.2p . we can express the time derivatives 
of $j as 

i$i,t = Qmi + 2iCQ$2 - Qx$2 
= - (Q*2), + 4iCQ$2 + 2Qi?$i 
= + 4iCg<f 2 + 2Qm,, (A.3) 

i$2,t = 2iCi?$i + - 4C'$2 - 
= + 2ic$2,. - 2i?g$2 

= $2,.x - 2($2,x - 2iC$2)$r'g'^2. (A.4) 

Because ^2 = 1^ and = Q', f !A.4p can be identified with the second 

equation in (]A.2|) . Combining ( lA.Sp with (11. lap and noting the relation 
^i,x = Q^2 (cf. f ll.2ap ). we can rewrite the time derivative of ^i^Q as 

= (<i>^i<i>i,,<i>ri)^ g + 2<i>r^$i,.$r^$i,.$r'g - ^i'q.. 

- 4ic$r'Q*2$r'g 
= - ($r'Q),, - 2$r'Q$2 (<f r'Q),. - 4ic$r'Q^2$r'Q- 

This verifies the first equation in (]A.2p . □ 

The correspondence relation flA.ip between the NLS system and the 
Chen-Lee-Liu system can be generalized for the higher /negative fiows of 
the hierarchies. In addition, it is easy to extend this idea to the case of 3 x 3 
(or even higher) matrix Lax representations, e.g., f l3.15p for the three-wave 
interaction system fl3.13p . 



B Semi-discrete Chen— Lee— Liu system 

In this appendix, we show that the Lax representation for the Ablowitz-Ladik 
lattice can provide the solutions of an integrable space discretization (semi- 
discretization, for short) of the Chen-Lee-Liu system. This result can be 
found by exploiting the transformation of dependent variables between the 
semi-discrete Gerdjikov-Ivanov system (14.1 7p and the semi-discrete Chen- 
Lee-Liu system. 
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Let Qn and Rn satisfy the nonreduced Ablowitz-Ladik lattice fl4.14p . Let 
and \l/2,n be the first and second components of an eigenfunction of the 
associated hnear problem fl4.15p . Then, we can state the following discrete 
analog of Proposition lA.ll 



Proposition B.l. The new pair of variables, 

-1 



satisfies a closed system, i.e., the semi- discrete Chen-Lee-Liu system [TH] 

(B.2) 

Moreover, multiplying \E'j_„. ^2/ factor e'^* m i/ie original definition of g„. a'^c? 
r„ in fIB.ll) . one can ac?c? t/ie terms +cqn and —crn to the left-hand side of 
equations flB.2p . The nonzero parameter C is nonessential, because it can be 
set as 1 by the transformation qn ='■ q'nC^^~^ ) ='■ '"nC"^" '^^^ redefinition 
of the parameters aC,"^ =: a' , bjQ =: b' . 



Remark. The index of the unit matrix J to indicate its size is suppressed 
in (1B.2P and below. For scalar g„ and r„, (lB.2p was previously studied 
in [S3HSS] (also see [SB])- Clearly, (]B.2p with a = — 6 = i is not an ideal 
semi-discretization of flA.2p . because it does not allow a local reduction to 
relate qn and r„ by complex/Hermitian conjugation. However, flB.2l) in the 
scalar case is an interesting system possessing an ultralocal Hamiltonian 
structure [SBHSS] , which plays a role in discussions on a tri-Hamiltonian struc- 
ture of the Ablowitz-Ladik lattice and related systems (cf. the concluding 
remarks in |79j). 

Proof. For brevity, we use the quantities $j^„ := j^nC"^^"^^^^^"^^^^ that 
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satisfy the following linear problem (cf. (I4.15P ): 



*2,n+l J 



1 p J_ r 



(B.3a) 



(aRn-l + ^Rn 



1 CaQn + f Q n-1 
—bRnQn-1 



(B.3b) 



we can rewrite the time derivatives of nQn-i and 

$2,n as 



Using flR3|l and fimj]) 

'^2,n as 



(C' - 1) a - (^^ - 1^ 6 ^I'nQn-l + a^^lQnRn-lQn-l 

b 

+ $1^.^ [aQn - bQn-2 - (a - - aQ^Rn-lQu-l + 6Qn-l-Rn-l<5n-2] 

= -c^a<i>i;^g„_i + a$i;ig„, (/ - c$2,n$i;^<5n-i) 

+ — $^^g„^i (/ - C$2,n$i;iQn-l) - &$i;nQn-2 + b^^^^Q n-1 Rn~lQn-2 

+ ^ (/ - C<f vign-l$2,n) '^li^Qn-l - b (^I - i$r,;,g„-l$2,„-l^ '^'M-lQn- 

+ 6 (/ - C$i:ign-l$2,n) (^<f riQn-l - $r,n-lQn-2) (B.4a) 
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-1 



and 

= a (C^$2,n - $2,n-l) (^I " ^$i;nQn-l$2,n-l^ 

+ b (^<|.2,„+1 - ^$2,n) (/ - C<^>r,nQn-l$2,n) , (B.4b) 

respectively. Because ^il^Qn-i = q-n and $2,ri = (IB.4p can be identified 
with (E2]). ' ' □ 

C (2 + l)-dimensional Chen— Lee— Liu systems 

In this appendix, we derive two (2 + l)-dimensional generahzations of the 
Chen-Lee-Liu system from the Lax representations for the two (2 + 1)- 
dimensional NLS systems considered in section [51 These results can be found 
by exploiting the transformation of dependent variables between each pair of 
the Gerdjikov-Ivanov system and the Chen-Lee-Liu system in 2 + 1 dimen- 
sions. 

First, we discuss the non-isospectral case considered in subsection 15. II Let 
Q and i?, together with the auxiliary fields / and g, satisfy the Calogero- 
Degasperis system (15. ip . Let and \E'2 be the first and second compo- 
nents of an eigenfunction of the associated linear problem (15. 2p . wherein the 
x-independent spectral parameter C, satisfies the non-isospectral condition 
C,t = 2CCy Then, we have the following proposition. 

Proposition C.l. The new pair of variables, 
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satisfies the (2 + 1)- dimensional Chen-Lee-Liu system 



m + KyQ + Qxy + 2iC^g + Aq^ - -q[B- {rq)y] = O, 

in + iCyV - r^y - 2iCrA + r^A + ^[B + {rq)y] r = O, 

A^ + [qr, A] = {qr)y, 
B^ = (-4iCrg + r^q - rq^)^ , 

, Ct = 2CCy, Cx = 0, 
where A and B are new auxiliary fields. 



(C.2) 



Remark. When C is a real function and the dependent variables are scalar, 
we can impose the complex conjugacy reduction r = iaq*, A* = —A, and 
B* = B with a real constant a. In particular, the reduction r = ig* simplifies 

dal to 

m + KyQ + qxy - 2Cqd-'dy{\q\') - 2qd-%{C\q\') 

+ iq^d-'dyHqf) - ^q [d^'dy {qlq - q*q^) - {\qW] = 0, (C.3) 



0. 



If C is a constant independent of y and t, ( IC.SP further reduces to (cf. (66) 
in [97] or (73) in [9H]) 

iqt + qxy-4^Cqd;;;^dy{\q\^) 

+ iq,d-'dyi\q\') - ^q [d-'dy (qlq - q*q.,) - i\qW] = 0. 



Considering a gauge transformation of the Lax representation [HS], we can 
relate this (2 + l)-dimensional Chen-Lee-Liu equation to the known (2 + 1)- 
dimensional Kaup-Newell equation [(4.3) in Ref . ISTj . 

Proof. For brevity, we use the quantities $j := \E'je^''^ {j = 1,2) instead of 
^j. From (l52iD, we obtain $1,^ = Q^2 and R = ($2,^ - 2iC$2) $r^- Thus, 
we can introduce the auxiliary field A as 

The compatibility condition of these "linear" PDEs for $1 implies the relation 

($r^g$2)^ -A.,- [<i>r'g$2, a] = o. (c.4) 
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Equation (IS.ldp for the auxiliary field g can be rewritten as 

g, = (-2iC$2$r'Q + ^2,.^i'Q)y ■ (C.5) 

Using f lS.lap and (I5.2bp . we can express the time derivatives of ^i^Q and $2 
as 

i ($r'Q), = -^I'Qxy + ^I'Qg - 2ic$r''^i,.'^r'Q + $r'Q.'^2$r'Q 



and 



-1 CI — ^ 1 ^ ' \ i ' y ^ ^ -v-^- 

2ic^$r'Q 



i$2,t = 2iC$2,3; + Ry^l " C/*2 

= 2iC$2,j; + [($2,x. - 2iC$2) -^T'I^ "^1 - 9^2 



(C.6) 



2,a;?/ 



2iC,<f2-($2,x-2iC$2)^-^?$2, 



(C.7) 



respectively. By setting $f ^g = q and $2 = ^1 we obtain a (2 + l)-dimensional 
system for q and r, but it appears asymmetric with respect to these variables. 
To rewrite it in a more symmetric form, we need only redefine the auxiliary 
field as ^ ^ 

g =: -iCyl + -{rq)y + -B. 

Thus, (IC.4p -f lU?7|l . together with the non-isospectral condition, verify (IC.2p . 
□ 



Second, we discuss the case without a spectral parameter, considered in 
subsection 15.21 Let Q and R, together with the auxiliary fields / and g, 
satisfy the Davey-Stewartson system (15. 5p . Let \E'i and \E'2 be the first and 
second components of a solution of the associated linear problem f l5.6p . Then, 
we have the following proposition. 



Proposition C.2. The new pair of variables, 



r 



*2, 



(C. 
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satisfies the (2 + 1) -dimensional Chen-Lee-Liu system 



' iqt + aq^x + hqyy + q \a{rq)^ + F] + 2bGqy = O, 
irt - ar^x - hvyy + [a{rq)^ - F]r + 2hryG = O, 



Fy + a {vyq - rqy)^ = O, 
^ + [qr, G] = {qr)y, 

where F and G are new auxiliary fields. 



(C.9) 



Remark. In the case of scalar variables, system (\C.9\i was studied in [901199] 
A special case of flC.9p for vector dependent variables was discussed in [100] 



Proof. From fl5.6al) . we have "ifi^x = Q^2 and R = \E'2,y\l'i ^. Thus, we can 
introduce the auxiliary field G as 

The compatibility condition of these "linear" PDEs for \E'i implies the relation 

(^r'Q^2)^ -G,- [^r'Q^2, G] = o. (c.io) 

Equation (15.5d|) for the auxiliary field g can be rewritten as 

gy = 2a{^2,y'^i'Q)^. (C.ll) 
Using (15 -Sap and ( ]5.6bp . we can express the time derivatives of "^i^Q and \I'2 



as 



i (^r'Q)i = -a'^i'Qxx - b^^'Qyy + ^^'Qg + b^^'^i^yy^^'Q 

+ a (vi/-igvi/2)2 v[/-ig - b (^r'Q),, - 2bG^^'Qy 
+ bG'^^^^Q + ^^'Qg + bG^^^'Q - a^r'Q^2,x-^r'Q 

- 2a^r'Q^2,x^r'Q + "^I'Qg - '2bG (^r'Q), (C.12) 
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and 



i^2,t = a^2,xx - - bmi^y + bRy^i 

= a^2,xx + b^2,yy " ^^2 - 2b^2,yG, (C.13) 

respectively. By setting ^l/^^Q = q and \l/2 = r, we obtain a (2 + l)-dimensional 
system for q and r, but it appears asymmetric with respect to these variables. 
To rewrite it in a more symmetric form, we redefine the auxiliary field as 

g =: a{rq)^ - F. 

Thus, flCnojl - flCnaD verify ^M)- □ 
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